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ABSTRACT 

Recent studies have shown that alpha Centauri B might be, from an observational point of 
view, an ideal candidate for the detection of an Earth-like planet in or near its habitable zone 
(0.5-0.9AU). We study here if such habitable planets can form, by numerically investigating 
the planet-formation stage which is probably the most sensitive to binarity effects: the mutual 
accretion of km-sized planetesimals. Using a state-of-the-art algorithm for computing the im- 
pact velocities within a test planetesimal population, we find that planetesimal growth is only 
possible, although marginally, in the innermost part of the HZ around 0.5AU. Beyond this 
point, the combination of secular perturbations by the binary companion and gas drag drive 
the mutual velocities beyond the erosion limit. Impact velocities might later decrease during 
the gas removal phase, but this probably happens too late for preventing most km-sized ob- 
jects to be removed by inward drift, thus preventing accretion from starting anew. A more 
promising hypothesis is that the binary formed in a crowded cluster, where it might have been 
wider in its initial stages, when planetary formation was ongoing. We explore this scenario 
and find that a starting separation roughly 15 AU wider, or an eccentricity 2.5 times lower than 
the present ones are required to have an accretion-friendly environment in the whole HZ. 
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1 INTRODUCTION 

The study of planet formation in binary stars is an important issue, 
as appro ximately ~ 20% of all de tected exoplanets inhabit such 
systems jPesidera & Barbieril2007h . Even if most of these binaries 
are very wide, at least 3 planets have been found in binaries of sepa- 
ration ~20 AU, the most interesting of them being the M > 1.6My„p 
planet at 2. 1 AU from the primary of the y Cephei sys tem, vyhich has 
been t he subject of several planet-formation studies ([Thebault et al.l 
|2004 iKlev & Nelsorill2008l : Ijang-Condell et al.ll2008l) ^ 

However, the system which has been most thoroughly investi- 
gated i s one where n o planet has been detected so far: alpha Cen- 
tauri Q.|Holman|&Wiegert ( 1997.) have shown that, for the copla- 
nar case, the < 3 A U region around a Ce n A can harbour planets on 
stable orbits, while iBarbieri et al] ( l2002h and lOuintana et al] ( l2002l 
|2007|) have found that the final stages of planetary formation, the 
ones leading from large planetary embryos to the final planets, are 
ossib le in the inner < 2.5 AU region. However, Thebault et al.l 
2008l hereafter TMS08) showed that the earlier stage leading, 
through runaway growth in a dynamically quiet environment (e.g. 
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iLissaueil 1993b , from kilometre-sized planetesimals to the embryos 
is much more affected by the companion star's perturbations. In- 
deed, the coupled effect of these secular perturbations and drag due 
to the gas in the nebulae leads to a strong orbital phasing according 
to planetesimal sizes. This differential phasing induces high colli- 
sion velocities for any impacts between objects of sizes si + S2, 
which can strongly slow down or even halt the accretion process. 
As a result, the region beyond 0.5 AU (0.75 AU for the most ex- 
treme cases explored) from the primary is hostile to mutual accre- 
tion of planetesimals, and thus to the in situ formation of planets. 
The presence of planets in these regions, however, cannot be ruled 
out, in the light of possible additional mechanisms not taken into 
account in the simulations, such as a change of the dynamical en- 
vironment once gas has been dispersed, or outward migration of 
planets formed closer to the star, or a change in the binary's sepa- 
ration during it's early history. 

We reinvestigate here these issues for the case of the secondary 
star a Centauri B. Indeed, a recent study by lOuedes et alj ( 1200 Sh 
has shown that, in addition to the final embryos-to-planets stage be- 
ing possible within ~ 2.5 AU from it, this star would be an almost 
perfect candidate for the detection of a potential terrestrial planet 
by the radial velocity method. This is mainly because of its great 
proximity, but also because it is exceptionally quiet. Moreover, be- 
cause of its lower luminosity, its habitable zone (HZ) lies much 
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closer, 0.5 to .9 AU jOuedes et al.l20o3) . than that of o- Cen A (~ 1 
to 1.3 AU, e.g. lBarbierietalj2002l) . 

In this letter, we first numerically study the dynamics of the 
planetesimal accretion phase in a similar way as for a Cen A, by es- 
imating the extent of the accretion-friendly region around the cen- 
tral star (section II). We then push our studies a step further than in 
TMS08, by quantitatively exploring two mechanisms which could 
potentially increase the odds for planet formation in the HZ (Sec- 
tion III). Conclusions and perspectives are given in the last Section. 



2 SIMULATIONS 
2.1 model and setup 

We consider a disc of 10"* test planetesimals extending from 0.3 to 
1.6 AU from the central star, starting on quasi-circular orbits such 
as initial encounter velocities (Ai'> < 1 m.s"', and with a physical 
radius j in the lkm< s < 10km ra nge. The binary's parameters are 
taken from lPourbais et al] ( l2002h : semi-major axis at = 23.4 AU, 
eccentricity = 0.52, Ma = IAMq and Mb = 0.93Mq. We fol- 
low the dynamical evolution of the planetesimal population for 
t = lO^'years, the typical timescale for runaway growth, under 
the coupled influence of the companion star's (here a Cen A) per- 
tur bations and gas drag, usi ng the deterministic code developed 
by iThebauIt & Brahid (|l9981. All test planetesimal collisions are 
tracked with the code's build-in close-encounter search routine, 
looking at each time step for all intersecting particle trajectories 
with the help of the classical "inflated radius" procedure (see sec- 
tion 2. 1 of TMS08 for more details). The gas disc is assumed ax- 
isymmetric with a volumic density profile following pj = pgo ''(\u)' 
where we consider as a nominal case the Minimum Mass Solar 
Nebula (MMS N), with q = -2.75 and p„o = 1.4 x IQ-'g.cm-^ 
( lHavashi|[T98lh . exploring pgo and q as free parameters in separate 
runs. 

All mutual encounters and corresponding encounter velocities 
are tracked and recorded in order to derive {Av)si,s2 matrices for 
all possible impactor pairs of sizes si and S2. These velocity values 
are then interpreted in terms of accreting or eroding encounters by 
comparing them, for each imp acting pair, to 2 threshold velocities 
(for more details, see lThebault et al.]|200dl2008l) : 

• ,„), the escape velocity of the impacting pair. This ve- 
locity gives approximately the limit below which single star-like 
runaway accretion is possible (when accretion rates are enhanced 
by the gravitational focusing factor). 

• the velocity above which impacts preferentially lead 
to mass loss rather than mass accretion. This value is obtained from 
estimates of the critical threshold energy for catastrophic fragmen- 
tation 2*(.ii.i2)- Given the large uncertainties regarding the possible 



values for Q* 



we consider 2 limiting values, v*, 



(s\ ,S2 )low 



and 



v*(.i],i2)high, corresponding to a weak and hard material assumption 
respectively. 

Finally, at any given location in the system, the global balance be- 
tween accreting and eroding impacts is estimated assuming a size 
distribution for the planetesimal population, weighting the contri- 
bution of each {su S2) impact by its normalized probability accord- 
ing to the assumed size distribution. As a nominal case, we assume 
a Maxwellian peaking at the median value s,„gj = 5 km, but other 
distributions (Gaussians, power laws) are also explored. 
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Figure 1. eccentricity vs. semi-major axis graph at the end (10''yeai's) of the 
nominal case run. The colour scale shows the different paiticle sizes. The 
two vertical lines indicate the estimated location of the habitable zone. 
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Figure 2. Relative importance of different types of collision outcomes, 
as a function of distance to the primary star (aCenB), at 1=10'' years for 
our nominal case. Each A\'(,| „) for an impact between bodies of sizes .vi 
and S2 is interpreted in terms of eroding impacts (if Ai' > v*(s| „)i,igi,), 
accreting impacts despite of increased velocities ("perturbed accretion" 
for iVjr(ji,j2) < < i'*(s|_„)io„), or unperturbed "normal" accretion 
(Av < vv.,r(.i,..i2))- For v*(sj,i2)io„ < Al' < v*(.,|_.,2)high, we do not draw any 
conclusions. The contribution of each impact is then weighted assuming 
that the size distribution for the planetesimals follows a Maxwellian cen- 
tered on 5km (see text for details). 



2.2 results: nominal case 

The (e,a) graph displayed in Fig[T] clearly illustrates the effect of 
differential orbital phasing as a function of sizes. Smaller planetes- 
imals align to orbits with lower eccentricities than that of the bigger 
objects. Ij. Note that in the outer regions (beyond ~ 1 AU), residual 
eccentricity oscillations are observed for the larger bodies. This is 
because these objects, which are less affected by gaseous friction. 



^ This differential alignment in e is strength ened by the associate d differ- 
ential alignment of the pericenter angle (e.g.. lThebault et al.l2006l) 
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Table 1. Outer limit racc{out) of the accretion-friendly inner zone, as defined 
by the region where the "green" and "blue" areas make up more than 50% 
of the collision outcomes (see Fig|2), for test runs exploring different free 
parameters: gas density pjo at lAU, slope cj of the pg cc gas density 
profile, and assumed size distribution for the planetesimal population 



Set-Up raccioM) 



Nominal case 


0.5 AU 


PgO = 1.4 X 10-"'g.cm-3 (O.lxMMSN) 


0.3 AU 


PgO = 1.4x 10-**g.cm-3 (lOxMMSN) 


0.8 AU 


Pg oc j-Q, viscous disc) 


0.45 AU 


Pg oc r 


0.4 AU 


Size dist.: Gaussian, cr^ = (2km)^ 


0.55 AU 


Size dist.: Gaussian, tr^ = (1 km)^ 


0.6 AU 


Size dist.: Gaussian, cr^ = (0.5 km)~ 


1.1 AU 


Size dist.: power law dN oc s^^ ^ds 


0.5 AU 



have not had the time yet to fully reach the equilibrium orbit forced 
by gas drag. 

Fig[2]shows how {Av>jj resulting from this differential phas- 
ing can be interpreted in terms of accreting versus eroding impacts. 
This graph displays the relative weights of the possible types of col- 
lision outcomes as a function of distance to o-CenB, for our nomi- 
nal case (MMSN gas disc), assuming that the planetesimal size dis- 
tribution follows a Maxwellian centered on 5km. By comparing this 
graph to the equivalent graph for a Cen A (Fig.3 of TMS08), we see 
that results are roughly comparable for both stars. Schematically, 
the region beyond ~0.5 AU from the primary is hostile to planetary 
accretion. This similarity of the results is mostly due to the fact that 
the mass ratio between the 2 stars is close to 1 (~ 0.85), and that 
the strongest perturbations of a Cen A on o- Cen B are partly com- 
pensated by the lower Keplerian velocities around that smaller star. 
One important difference, however, is t hat the HZ around a Cen B 
lies closer to the star: 0.5 to 0.9 AU jOuedes et al]|200 8'). This 
means that its innermost part lies at the outer-edge of the accretion- 
friendly region, thus opening a possibility, although marginal, for 
planet formation there. 

2.3 parameter exploration 

To check the robustness of these results with respect to our choice 
of parameters for the nominal case, we firstly explore other possible 
gas disc profiles, by varying both the slope q and the density pgo at 
lAU. We find that the only case which leads to a wider accretion- 
friendly zone is that of a dense lOxMMSN gas disc, for which this 
region extends up to ~ 0.8 AU (Tab[T]l. All other gas-disc cases 
lead to equally narrow or even narrower accretion-friendly regions 
than for our nominal case, except for a purely academic gas-free 
case. Of course, our axisymmetric gas disc assumption is probably 
a crude simplification of the real behaviour of a gas disc perturbed 
by a compani on star. However, prelimi nary Hydro-l-N-body runs 
performed by Paardekooper et all j2008l) show that the build up of 
an eccentricity for the gas-disc, altough this mechanism is yet not 
fully understood, always lead to higher perturbations and impact 
velocities in the planetesimal population. As a consequence, our 
axisymmetric disc case might be regarded as a best-case scenario 
and our results to be on the conservative side (for a more detailed 
discussion and justification for our choice of an axisymmetric disc, 
see TMS08). 

Another crucial parameter which has been explored is the as- 
sumed size distribution for the planetesimal population, as it is a 
parameter poorly constrained from planet-formation models (see 



discussion in iThebault et alj|2006l) . For any "reasonable" distri- 
butions, we do not observe a significant displacement of the ac- 
cretion/erosion radial limit. Only a very peaked, almost Dirac-like 
Gaussian of variance cr^ = (0.5 km)^ leads to mostly accreting im- 
pacts in the whole HZ (Tab[l). 



3 DISCUSSION 

The results of the previous section show that, as was the case 
for a Cen A, the region allowing km-sized planetesimal accre- 
tion is much more limited than the one allowing the final stages 
of planet formation (starting from large embryos): the accretion- 
friendly zone only extends up to ~ 0.5 AU for our nominal case, as 
compared to the esti mated 2.5 AU for the embryos-to-planets phase 
jOuedes et al.ll2008l . and reference therein). 

The only cases leading to an accretion-friendly environment 
in a wider area than for the nominal run are: a) the gas free case, b) 
the very-peaked Gaussian size distribution, and c) the lOxMMSN 
run. Cases a) and b) are obviously very unrealistic test runs, but 
even case c) is probably here of limited physical significance. In- 
deed, assuming such a high-mass disc around a 0.9Mo star appears 
as a rather extreme hypothesis, whic h is not backed by theoretica l 
estimates or observational data (e.g. [Andrews & Williamsll2007h . 
For a binary system, this issue gets even more critical, as the ex- 
pected outer truncation of the circumprimary gas disc probably 
de pletes it from a subst ancial amount of it mass reservoir (e.g., 
Ijang-CondeU et alj|2008l) . making the lOxMMSN hypothesis even 
more questionable, although it cannot be completely ruled out. 

Apart from this massive gas disc case, the r„„(oBr) ~ 0.5 AU 
limit allows in principle for planetesimal accretion to occur in the 
innermost part of the habitable zone, which is estimated to extend 
from 0.5 to 0.9AU (as compared to ~ 1 to 1.3 AU for o-CenA). 
However, even if some planetesimal accretion is possible, it could 
probably not be what it is in the standard planet formation sce- 
nario. Indeed, even if the ~ 0.5 AU region is accretion-friendly, 
relative velocities there have values higher than what they should 
be in an unperturbed system (t hat is, approxim ately the escape ve- 
locities of km-sized bodies, e.e l Lissaueill99"3h . This means that the 
gravitational focusing factor, which is the cause of the fast runaway 
growth mode, is significantly reduced or cancelled. Growth would 
thus have to be m uch slower, maybe orderly or of the "type 11" run- 
away identified bv lKortenkamp et^ ( l200ll) . In short, planetesimal 
accretion is not possible in the HZ, except in the innermost region 
around 0.5 AU where it should be significantly different from what 
it is around a single star. 

However, no conclusions regarding the presence of planets in 
the HZ can be directly drawn from these results, because of the 
possible effect of additional processes not taken into account in our 
simulations. In TMS08, we reviewed some of these possible mech- 
anisms, and concluded that, while some of them might not signifi- 
cantly affect our conclusions (like bigger "initial" planetesimals or 
re-accretion in the subsequent gas-free disc), others might indeed 
open for a possible presence of a planet in the terrestrial region. 
We reinvestigate here these mechanisms in more details, focusing 
specifically on the two most promising scenarios, that of orbital 
re-phasing during an extended gas dispersion phase and that of a 
greater separation for the early binary. 

3.1 re-phasing during gas dispersion 

Ixie & Zhou! ( |2008|) have investigated, for the specific case of the y 
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Figure 3. (e,a) graph, after ; = 2 X 10 years, for the run with gas dissipation. 
Gas dissipation starts at t=10'*years and has a timescale Tf/,,, = lO^years 
(see text for details) 



Figure 4. Dominant collision outcome (erosion, perturbed accretion, un- 
perturbed accretion, see Fig|2j in the habitable zone (0.5-0.9 AU) around a 
Centauri B, for different configurations of the binary's orbit. 



Cephei binary, what happens to a planetesimal swarm during the 
gas dispersal period. They have identified an interesting mecha- 
nism: if the gaseous component is slowly removed from the disc, 
all planetesimal orbits are progressively re-phased towards the 
same orbits regardless of their sizes. This leads to lower, possi- 
bly accretion-friendly impact velocities. They concluded that this 
might allow, under certain circumstances, for planetesimal accre- 
tion to start anew once this realignment is achieved. We explore this 
possibility for the present o rCenB case, assumi ng the same param- 
eters for gas dispersion as IXie & Zhoul ( l2008h : Pg {t It diss) , 
with Trfijj = 10''years. The dissipation is started at lO^'years, the 
end of our nominal run, when all orbits have reached their size- 
dependent alignment in the HZ. 

Fig[3] shows the dynamical state of the system after 2x^,5, = 
2xl0'years and clearly illustrates the efficient re-phasing of all or- 
bits, which is almost complete in the > 0.8 AU region. In this outer 
region, (Av) are low enough to allow accretion of all remaining ob- 
jects. However, this encouraging result is undermined by several 
problems. The first one is that, in this accretion-friendly outer re- 
gion, no object smaller than ~ 4 km is left after 2x^,5,, because the 
time it takes for most of the gas to be dispersed is long enough 
(a few 10' years) to have a significant inward drift of all smaller 
bodies!^. True, in most of the habitable zone small planetesimals 
are still present, but here full orbital re-phasing is not achieved yet, 
and the dynamical environment, even if the situation is slighly im- 
proved with respect to the nominal run, is still globally hostile to 
accretion. Only after ~ STjiss = 5 x lO'years do (Ai'> become low 
enough to allow accretion in the whole HZ, by which time all small 
objects have here also been removed. 

The second problem is that, even for the larger > 4 km 
planetesimals, the favourable conditions shown in Fig[3]are only 
reached after an accretion-hostile transition period of a few r^,,,. 
The question is then how these objects might survive this long pe- 
riod during which most impacts will erode them. These erosion 
processes cannot be followed with the N-body models used here, 

' This issue had already been identified bv lxie & ZhoiJ jlOOSi) but didn't 
show up in their simulations because of the numerically imposed re- 
injection at the outer boundary of all objects lost at the inner one 



but it is likely that most large planetesimals will be fragmented into 
smaller debris which will be quickly removed from the system by 
fast inward drift. 

For all these reasons, we do not believe that gas dispersal 
opens the possiblity for planetesimal accretion to start anew after 
a few rj,,j. How ever, the re-p hasing mechanism during gas disper- 
sal identified bv lXie & Zhoul 12008) is definitely worth studying in 
more details. One crucial, and yet almost unexplored issue is for in- 
stance the relative timing between the planetesimal accretion phase 
(in particular its start) and the parallel viscous evo lution of the gas 
disc. Indeed, as pointed out bv lXie & ZhoJ j2008l) , the steady evo- 
lution of the gas disc profile, due to viscous angular momentum 
redistribution and accretion, can lead to a gas density decrease 
of several orders of magnitudes long before its "definitive" and 
abrupt dispersion, at 5 to lOMyrs, resulting from photoevaporation - 
induced disc truncature (see e.g.. Fig. 1 of [Alexander et al]|2006h . 
The question is where (or when) to place the planetesimal accre- 
tion phase in this picture. 

3.2 Wider initial binary separation 

Most stars are bom in clusters. As such they spend their early life 
in an environment where the risk for close stellar encounters is rel- 
atively high. For a binary system, such encounters could have dra- 
matic effects, breaking it up or significantly modifying its initial 
orbit. Conversely, a present day binary cou ld be what is lef t of an 
initially triple (or more) system (Marzari & Barbierill2007al lbl). As 
a consequence, there might be an important difference between the 
present day orbit of a binary system and the one it had when plan- 
etary formation processes were at play. 

In order to have an idea of the statistical chances for an a 
Centauri-like binary to have an initially diff^eren t orbit, we take as 
a reference the work of iMalmberg et al] ( |2007|) . who studied the 
dynamical evolution of a typical open cluster. They found that most 
binaries of initial separation < 200 AU are not broken up during 
the cluster's lifetime. However, a fraction of these "hard" binaries 
suffer stellar encounters which alter their orbits. The net result of 
these encounters is to shrink the binaries' orbits, so that separations 
of young binaries tend to be larger, with an envir onment which is 
thus more favourable to planet formation. Fig.4 of lMalmberg et al.] 
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l l2007h indicates tliat ~ 50% of binaries with separation ~ 20AU 
have suffered an orbit-modifying encounter. There is thus roughly 
1 chance in 2 that a Centauri had a wider separation during the 
early stages of planet formation if it formed in a dense cluster. 

The crucial question is here which initial orbital configura- 
tions could have allowed planetesimal accretion to occur in the hab- 
itable zone of a Cen B. We investigate this issue by running a series 
of 31 test runs for 20 < Ofc < 50 AU and < < 0.8 (FigH. Not 
surprisingly, we see that present day a Cen B lies almost at the limit 
where accretion, although highly perturbed, becomes possible in 
the innermost part of the HZ. For the same ei, = 0.52 as today, par- 
tial perturbed accretion in the HZ is possible for 23 < < 37 AU. 
In order to have the HZ become fully accretion-friendly, one needs 
fli, > 37 AU, while unperturbed single star-like accretion requires 
flj, > 52 AU, more than twice the present day separation. The sit- 
uation gets of course more favourable if one assumes a lower ini- 
tial et for the binary. As an example, with et = 0.26 (half of the 
present day value), then the whole HZ becomes accretion friendly 
for ai, > 26 AU and allows unperturbed accretion for at, > 40 AU. 

These results are of course only a first step towards under- 
standing this complex issue. A crucial point is to quantitatively es- 
timate what the probability for a given Aaj and Aej "jump" are. As 
an example, how likely is it for a present day a Centauri-like system 
to have originated from the accretion-friendly regions (3) or (4) in 
FigEf The answer to these questions requires very detailed and en- 
compassing numerical explorations which have, to our knowledge, 
not been performed yet. Marzari & B arbieri ( 2007a, b) have studied 
in great detail the dynamical outcomes of individual encounters, 
but decoupled fro m the general cluster c ontext, while the cluster 
evolution study of iMalmberg et al.] ( l2007h does not give statistical 
information about the amplitudes of binary orbital changes. 

4 CONCLUSIONS AND PERSPECTIVES 

We have numerically investigated planetesimal accretion around a 
Centauri B. Our main conclusions can be summarized as follows: 

• Planetesimal accretion is marginally possible in the innermost 
parts, ~0.5 AU, of the estimated habitable zone. Beyond this point, 
high collision velocities, induced by the coupling between gas fric- 
tion and secular perturbations, lead to destructive impacts. More- 
over, even in the ~ 0.5 AU region, (Av) are increased compared 
to an unperturbed case. Thus, "classical", single-star like runaway 
accretion seems to be ruled out. 

• These results are relatively robust with respect to the planetes- 
imal size distribution or gas disc profile, except for a very massive, 
and probably unrealistic lOxMMSN gas disc. 

• We confirm the conclusions of several previous studies that 
the planetesimal-to-embryo stage is much more affected by binarity 
effects tha n the subsequ e nt em bryos-to-final-planets stage. 

• As in IXie & Zhoul ( l2008l) . we find that later progressive gas 
dispersal reduces all (Av> to values that might allow accreting im- 
pacts. However, we find that the system has first to undergo a long 
accretion-hostile transition period during which most of the smaller 
planetesimals are removed by inward drift and most bigger objects 
are probably fragmented into small debris. Thus, the positive effect 
on planetesimal accretion is probably limited. 

• We quantitatively investigate to what extent a wider initial bi- 
nary separation, later reduced by stellar encounters in an early clus- 
ter environment, could have favoured planetesimal accretion. We 
find that, for a constant e^, an entirely accretion-friendly HZ re- 
quires an initial at, > 37 AU, while normal unperturbed accretion 



is only possible for ai, > 52 AU. The statistical likehood of such 
orbital changes in early open clusters remains to be quantitatively 
estimated. 

We conclude that, although the presence of planets further 
out cannot be fully excluded, the most likely place to look for an 
habitable terrestrial planet would be around 0.5 AU. According to 
lOuedes et alj ( |2008[) . the detection of such a planet could be pos- 
sible, if its mass is > l.SM^, after 3 years of high cadence obser- 
vations, provided that the noise spectrum is white. Let us however 
point out again that the ~ 0.5 AU region around o-CenB is dy- 
namically perturbed by the binary, so that the formation of a planet 
there cannot have followed the standard single-star scenario (unless 
the initial binary separation was much larger, see Section 3.2). The 
planet formation process in such a perturbed environment, which 
could concern many binary systems of intermediate (~ 20 AU) sep- 
arations, is an important issue which remains to be investigated. 
Another crucial issue which remains to be investigated is what hap- 
pens before the phase studied here. Indeed, our study implicitly as- 
sumes that km-sized planetesimals could form from smaller grains 
and pebbles, but the validity of this assumption has to be critically 
examined. 
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